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Abstract. The hydrogen bond is a fundamental ingredient to stabilize the DNA and RNA 
macromolecules. The main contribution of this work is to describe quantitatively this interaction 
as a consequence of the quantum confinement of the hydrogen. The results for the free and 
confined systems are compared with experimental data. The formalism to compute the energy 
gap of the vibration motion used to identify the spectrum lines is the Variational Method allied 
to Supersymmetric Quantum Mechanics. 


1. Introduction 

One of the most important interactions to stabilize the DNA and RNA macromolecules is the 
hydrogen bond between pairs of nitrogenous bases. It promotes a link among the complementary 
pairs of bases, maintaining the double helix structure. In the DNA the pairs of bases are 
guanine-cytosine and adenine-thymine; in the RNA the thymine is substituted by uracil, [Ij. In 
a recent work, the hydrogen bonds present in proteins and in smaller molecules were studied 
by considering the hydrogen quantum confinement in the NH and OH molecular groups, [2], In 
the present work, we extended the results for the nucleic acids, which present the NH molecular 
group. 

The hydrogen bond is present in various systems, [3], and it is classified as a donor-acceptor 
interaction. From the experimental point of view, it is possible to identify if a particular 
molecular group is forming an hydrogen bond due to the change in the vibrational spectrum. 
This change can be observed by different experimental techniques, as for instance the infrared 
spectrum analysis. 

The main point of this work is to show that the hydrogen bond can be treated as the quantum 
confinement of the hydrogen, which would change the transitions related to the NH group present 
in nucleic acids. As usual, the Morse potential, pt], is used to simulate the vibration of the NH 
molecular group, m-m- The parameters used are the ones found in literature for these groups 
without the hydrogen interaction, i.e., without confinement. The geometrical restrictions due to 
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Figure 1. Schematic one-dimensional model, where x c is the covalent radius of the atom 
donor, x vw is the van der Waals radius of the atom and acceptor and x eq is the equilibrium 
distance of covalent bond. X is the donor atom and Y represents the atom acceptor. 


the hydrogen bonds are obtained from standard experimental results for this kind of systems, 
[3j. The calculation of the energy levels is made through the variational method with test 
functions obtained from Supersymmetric Quantum Mechanics formalism, SQM. This procedure 
has alrealdy been used in different quantum problems, with and without confinement, m-m, 
including specific studies involving the Morse potential, m , which corroborates its applicability 
in the present case. 

In section 2 the adopted model for the hydrogen bonds is presented. We then present a review 
of the SQM formalism in section 3. In Section 4 the specific mathematical formalism concerning 
the Morse potential is developed, including the identification of the wave functions to be used 
in the Variational Method. The numerical results, with the appropriate parameters of the NH 
group are presented as well as the comparison with the experimental data are in section 5. In 
section 6 the conclusions are discussed. 

2. Model 

Figure 1 shows a schematic model of a hydrogen bonding between the donor atom X and the 
acceptor atom Y, here denoted by X-H- • -Y. The hydrogen bond is represented by H- • -Y and 
covalent bond is represented by X-H. The X and Y atoms are considered as hard spheres and the 
hydrogen is treated as a material point, since the electron cloud is shared with the donor and the 
acceptor atoms. The origin of the coordinate system is located at the center of the atom X. Any 
distance used is computed from this point and it is considered center-to-center. The covalent 
radius is adopted by the condition of impenetrability of the donor atom X; the employment 
of the van der Waals radius is justified because it is a measure of the volume excluded by the 
acceptor atom Y, i.e., the hydrogen cannot penetrate this region. 

Thus, the hydrogen nuclei may oscillate within a region determined between a minimum 
distance ( x m .i n ), given by the covalent radius (x c ) of atom X, and a maximum distance (x max ), 
given by the length of the hydrogen bond ( l x — y ) subtracted by the van der Waals radius (x vw ) 
of the atom Y, i.e.: 

Xmax = lx-Y X vw ( 1 ) 

The vibrational energy of this quantum system is suitably described by the confined one¬ 
dimensional Morse potential, widely known to well describe oscillations of diatomic molecules, 
m, m-m, m- The energy absorbed/emitted is obtained from the difference between the 
energies of the ground state and the first excited state. Since the confined Morse potential 







is non-exactly solvable, the energies of those states evaluated by using the variational method 
associated to the formalism of SQM. 

3. Supersymmetric Quantum Mechanics, SQM 

SQM is a well known algebraic method commonly used to exploit different aspects of non- 
relativistic quantum mechanical systems. It is particularly efficient to solve exactly all the 
quantum potential problems by means of the Hamiltonian factorization, thus providing the 
entire spectrum of wave functions and respective energies and a hierarchy of Hamiltonians, all 
related by the supersymmetric algebra, [Ill- 

Using the superalgebra, a given Hamiltonian H\, which has the energy E q 1 "* and wave function 

1 I'q 1) (x) for the lowest state, can be factorized in terms of bosonic operators Af, such that the 
Schrodinger equation can be rewritten (in h = c = 1 units) as 

Hi = ~£i + M*) = AUl + E o ] (2) 

where 

Hr^ 1) (x) = 0. (3) 

The following realization of the bosonic operators in terms of the superpotential w\(x) 

if= (4 + ' ,iW ) (4) 


results in the Riccati equation 


w l - ~^ w i = Vi(x) - ( 5 ) 

and defines the ground state wave function in terms of the superpotential 

Tg 1 ^) = Nexp(— [ wi(x)dx). (6) 

Jo 

It can be shown that a n-members hierarchy of Hamiltonians can be generated, all related 
by the superalgebra, m, 

H n = A+A~+E^ , ( 7 ) 

where 

A n = (t^ + w„(x)) (8) 

with simple relations connecting the eigenvalues and eigenfunctions of the n-members, 

ti 11 = AtA}...^ , EW = 4" +1 >. (9) 

In equation Q, the upper index between parentheses refers to the hierarchy Hamiltonian 
member and the lower index to the level within the hierarchy. 

At this point, it should be remarked that if the potential is non-exactly solvable, such as in a 
confined system, an approximation technique is needed and the variational method has already 
appeared as fully appropriate, [12]. 


4. Mathematical Formalism 

The one-dimensional Morse potential, [3], defined in terms of the parameters D e , the dissociation 
energy of molecule; /?, the parameter related to the width of potential well and x eq , the inter- 
nuclear equilibrium distance, is given by 

V M = D e {e~ 2 ^ x ~ x ^ - 2e~^ x ~ XEq) ). (10) 


By treating the molecule as an oscillator of frequency v, we can Taylor expand the Morse 
potential around the equilibrium position to find, in second order, the /3 parameter given by 


P = ( 
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( 11 ) 


and fix it through experimental data. With the change of variables, y = (3x, the dimensionless 
Schrodinger equation for a molecule subject to this potential is given by 


with 


and 
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where y is the reduced mass of the system. The Hamiltonian operator in equation (12) is given 

by 2 

Hm = —pr + A 2 ( e - 2 ^-^) - 2e~^ y ~ Veq> ) (15) 

dy 1 


Since the Schrodinger equation is exactly solvable, we identify Hm with H\ of equation Q and 
proceed the process of factorization, calculating the whole hierarchy. The result is, |12| 


V n+ i (y) = A 2 (e 2(y Veq) -2e ^ Veq) ) + 2nAe ^ Veq) 


w n+ i(y) = -Ae-<»-»"> + (A - 

e („ + l ) = _( A _2n + l )2i (16) 

In addition we can calculate exactly the wave functions. The ground state wave function is 
evaluated by using equation ([6]) and is given by 

^ ( o\y) oc exp ^—Ae~^ y ~ yeqS> — y (A — -)^ (17) 

This wave function is then used to evaluate, through the SQM algebra, the wave function for 
the excited states by using equation (9); the first excited state is given by 

T^(y) oc —2 — A + l^j exp (—Ae^^ y ~ yeq ^ — y(A — —) 


(18) 









Figure 2. Schematic model hydrogen bonds present in the AT pair in DNA. 


The wave functions given in © and ( fT8] ) are the ground state and the first excited state 
wave functions of the exact system, non-confined or with no hydrogen bond. We refer to 
them as T n(non-c<mfined)i n = 0,1. As our aim is the lowest state and the first excited 
state of the confined system, with the hydrogen bond, a potential barrier is introduced, i.e., 
Vi {y m ,ax) = hi (Umin) = oo, meaning that these wave functions must be zero outside the barrier. 
Consequently, the probability of finding the hydrogen outside the confinement region must be 
zero. 

As the system is no longer exact, it is necessary to use an approximative method, which 
we choose to be the variational method associated to SQM. From the SQM point of view, 
this is equivalent to introduce two additional terms in the superpotential w\, associated to the 
superalgebra, [T3] , [13] . Adopting this procedure, the wave functions to be used in the variational 
method are given by 


^n\confined)°^^n\non-confined)-(y™ax J/)G/ Umin ) Tl — 0,1 , 
with the corresponding energies obtained by solving 


E n — 


iy::: khm^v 
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(19) 


( 20 ) 


where Hm is given by equation (15) and are the trial wave functions given by equation (19), 
which depend on © and ( |18[ ). The integration limits are the covalent radius of the atom for 
which the hydrogen is chemically bound ( x m i n ) and the van der Waals radius complementary of 
hydrogen bond (x max )- The results of energy are presented in cm -1 units (inverse wavelength), 
the same as the ones from spectroscopic results. 


5. Results 

The hydrogen bonds studied involve the NH group. The necessary parameters, presented in 
Tables 1 and 2 below, were obtained experimentally, as indicated. 

Table 1. Parameters used for the calculations of energy to the NH group, m- 


NH 

Equilibrium distance (x eq ) 

1.0362 A 

Harmonic vibrational wavenumber (a;) 

3282cm -1 

Energy dissociation ( D e ) 

3.47 eV 

Reduced mass (fi) 

1.5614.10 ~ 27 kg 
















Figure 3. Schematic model hydrogen bonds present in the AU pair in RNA. 


Table 2. Covalent radius and van der Waals radius used for the nitrogen and oxygen atoms 

m , M- 


_ N O 

Covalent radius (x c ) 0.71 A 0.66 A 

van der Waals radius (x vw ) 1.55A 1.52 A 


5.1. Results for NH group 


5.1.1. Non-confined case 

Using the values given in Table 1 the value of the Morse potential (5 parameter, determined 
by the relation ( |ll| ) for NH group is /3 = 2.3168.10 10 m -1 . The value of the constant A for NH 
group, given by equation (13), is A = 17.0534. 


In the situation free of any confinement the calculation of the transition from the first excited 
state to the ground state is made according to equation (16). In this case, in spectroscopic units, 
this value corresponds to E\ — Eq = 3089cm -1 and should be found in spectroscopy. In studies 
by infrared spectroscopy with polypeptides it is observed an absorption band in 3090cm -1 and 
it has been described as a vibration frequency of the NH molecular group, |21j . This value is 


found very close to the calculated and can be identified, according to the model proposed, as 
vibration of NH without the formation of hydrogen bonding. 


5.1.2. Confined case 

When the NH group form molecular hydrogen bonding the vibration frequency is shifted. An 
interesting system where you can observe this effect are in the nucleic acids, DNA and RNA, 
because between the base pairs there are hydrogen bonds that allow the maintenance of the 
structure of the macromolecule. Here we analyze the vibrations of the NH group in hydrogen 
bonding between base pairs adenine-thymine (AT) and adenine-uracil (AU). In these base pairs 
two hydrogen bonds are formed, one with oxygen and the other with nitrogen. The schematic 
models are shown in Figures 2 and 3. 

The lengths of the hydrogen bonds of AT and AU present in the DNA and RNA molecules, 
respectively, are well known in the literature and are obtained by X-ray diffraction, [22]. The 
distances are given between the donor atom and acceptor atom, because by X-ray diffraction 
it is not possible to determine the hydrogen positions. The lengths of the hydrogen bonds are 
shown in Table 3 below. 







Table 3. Length of hydrogen bonds in nucleic acids, [ 22] . 



AT 

AU 

(N...N 

2.927 A 

2.820A 

(N...O) 

2.861 A 

2.950A 


These lengths are used in equation ([Tj) to determine the hydrogen maximum limit of 
oscillation. 

Solving equation (20) through the variational method, we found the transition energies for 
each bond; they are summarized in Table 4. These values are compared with results given in 
literature obtained experimentally. 


Table 4. Energy levels of the NH group in the hydrogen bonds between the base pairs 
AT and AU. Eq is the ground state energy, E\ is the energy of the first excited state and 
AEi^o = E\ — Eq. in (cm' 1 ) units. 





E 0 

Ei 

A-Ei-s.0 

AT 

(N- ■ 

■■N) 

-26330 

-23122 

3218 


(N-. 

•o) 

-26314 

-23039 

3275 

AU 

(N- ■ 

■■N) 

-26225 

-22794 

3431 


(N- ■ 

•o) 

-26340 

-23175 

3165 


Femtosecond vibrational spectroscopy with DNA oligomers realized with an alternating 
sequence of 23 adenine and thymine base pairs, [23], shows that there are three absortion peaks 
in the vibrational spectrum interesting to remark. 

One is between 3000 and 3100cm -1 and two around 3250cm -1 , one being located between 
3200 and 3250cm -1 , related to the N...N hydrogen bond in the AT pair, and the other between 
3250 e 3300cm -1 , related to the N...0 hydrogen bond in the AT pair. These last two peaks 
can be identified, according to our model, as the vibration of the hydrogen bond in NH, with 
a deviation between 0.56%-0.98% and 0.76%-0.77% respectively to the two ranges. The first 
peak, located around 3000 and 3100cm -1 , was found to be in fact at 3075cm -1 , and can be 
associated, according to our model (3089cm -1 ), with NH group without hydrogen bond, with a 
deviation of 0.46%. 

Thus, these experimental results of the AT of DNA and AU of RNA oligomers show perfect 
agreement to the results obtained through our method. 

Other infrared spectroscopy studies with RNA oligomers, |24j . show the vibration spectrum 
when the NH participates in hydrogen bonds with some explicit numerical values. One of the 
peaks of the spectrum is located at 3185cm -1 , very close to 3165cm -1 relative to the binding 
(N- • -O) calculated by our method, with a deviation of only 0.63%. A peak between 3050 e 
3100cm -1 can also be observed, which can de associated to the non-confined situation, with no 
hydrogen bond of NH group of our model, (3089cm -1 ). 

6. Conclusion 

In this work the quantum confinement of hydrogen present when hydrogen bonds are stablished 
was studied through the variational method associated to SQM. The vibration spectrum of NH 
group was evaluated in the situations of non-confinement, i.e., with no hydrogen bonds and 
when there is confinement, i.e., when the hydrogen bonds are formed. We focused our attention 
to the case of DNA and RNA oligomers, which present the NH group. As a result we conclude 
that the presence of hydrogen bonds promote significant changes in the vibration spectrum of 
NH group. The comparison with experimental data reported in the literature showed perfect 











agreement, with percentage errors smaller than 1%, which reenforces the applicability of the 
proposal. 

We conclude that the variational method associated to SQM was successfully applied to the 
calculation of the vibrational spectrum of NH molecular group with the principle of quantum 
confinement of hydrogen bonds. The treatment proposed here can be extended to the study of 
other molecular groups, as FH and CH, involved in hydrogen bonds. 
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